Four-dimensional CP^+ U(l) lattice gauge theory for 3D antiferromagnets: 
Phase structure, gauge bosons and spin liquid 
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In this paper we study the lattice CP^ model in (3 + 1) dimensions coupled with a dynamical 
compact U(l) gauge field. This model is an effective field theory of the s — ^ antiferromagnetic 
Heisenberg spin model in three spatial dimensions. By means of Monte Carlo simulations, we 
investigate its phase structure. There exist the Higgs, Coulomb and confinement phases, and the 
parameter regions of these phases are clarified. We also measure magnetization of 0(3) spins, 
energy gap of spin excitations, and mass of gauge boson. Then we discuss the relationship between 
these three phases and magnetic properties of the high-Tc cuprates, in particular the possibility of 
deconfined-spinon phase. Effect of dimer-like spin exchange coupling and ring-exchange coupling is 
also studied. 
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The CP^ spin model plays an important role in vari- 
ous fields of physics not only as a tractable field-theory 
model that has interesting phase structure, but also as 
an effective field theory for certain systems in condensed 
matter physics and beyond. In particular, the CP^ model 
corresponds to the Schwinger-boson representation of the 
s = \ antiferromagnetic (AF) quantum spin model, i.e., 
the AF Heisenberg model[l]. The CP^ model is much 
more tractable than the original AF Heisenberg model, 
and its phase structure and critical behavior have been 
investigated both analytically and numerically. The sys- 
tem intrinsically contains compact U(l) gauge degrees of 
freedom, and their dynamics determines the low-energy 
excitations in AF magnet. That is, if the gauge dynamics 
is in the deconfined- Coulomb phase, the low-energy ex- 
citations are the = | spinous. On the other hand, the 
Higgs phase corresponds to the Neel state with a long- 
range AF order, and the confinement phase is a valence- 
bond solid (VBS) state in which spin-triplet low-energy 
excitations appear. 

Most of the previous studies exploring a possible 
deconfined-spin-liquid phase have considered the two- 
dimensional (2D) (doped) AF Heisenberg model or its 
path-integral representation, the three-dimensional (3D) 
CP-^ model. In these cases, the Coulomb phase may be 
possible if there exist sufficient number of gapless matter 
fields that couple to the gauge field. In fact, the problem 
whether deconfinement in a 3D U(l) gauge theory is pos- 
sible was once a matter of controversy, but it seems now 
settled down, i.e., when the density of gapless excitations 
in matter sector is sufficiently large, the deconfinement 
phase takes place as a result of shielding mechanism [2]. 

In the present paper, we shall consider a higher- 
dimensional version of the 3D CP^ model, i.e., the 4D 
CP^ model coupled with a dynamical U(l) gauge field. 
This 4D CP^+U(1) gauge model is viewed as an effective 
field theory of the 3D AF Heisenberg model. From the 
gauge-theoretical point of view, the deconfinement na- 



ture is enhanced in (3 + 1) D case because the Coulomb 
phase exists even in the pure 4D U(l) gauge system that 
involves no matter fields in contrast to the pure 3D U(l) 
gauge system that has only confinement phase. There- 
fore, it is interesting to study the phase structure of this 
4D CP^ gauge model. We shall first consider the CP^ 
model for the 3D AF Heisenberg model with uniform 
nearest-neighbor spin coupling and then the CP^ model 
for the 3D AF Heisenberg model with nonuniform dimer- 
like coupling and ring-exchange spin coupling. 

The CP^ variable Zx is put on the site x of the 4D 
hypercubic (space-imaginary time) lattice. Zx is a two- 
component complex field satisfying the CP^ constraint, 
Zx = {zxi,Zx2)\ Z)a=i \^xa\^ = ^' path-iutcgral, the 
s = ^ spin operator Sr at the spatial site r is mapped to 

a classical 0(3) spin Sx = ZxOZx satisfying Sx • Sx = ^ 
where x = (r, xq) and a are the Pauli spin matrices. The 
U(l) gauge field Ux^ = exp(z6>^^) (/i = 1, • • • , 4 is the 
direction index and denotes also the unit vector in the 
/i-th direction) is put on the link (x, /i) connecting sites 
X and X -\- jii. 

The CP^ field theory in the continuum is derived from 
the 3D AF Heisenberg model by integrating out the half 
of the CP^ variables (on all the odd sites) by assuming 
a short-range AF order [1]. In order to study the model 
numerically, we reformulate the model by putting it on 
the 4D hypercubic lattice, which has the correct contin- 
uum limit. Then the action S of the CP^ model is given 



S — — ^ ^ {^Zx-\-^^aUx^Zxa + H.I 

{pxvUx+v.^Ux+^i^vUx^ ^^'^^ (1) 
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where ci and C2 are parameters of the model. 

Qualitative estimation of the parameters in the action 
^ is obtained as follows. By choosing the lattice spacing 
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FIG. 1: Phase diagram of the 4D CP^ model {[T]) in the C2 -ci 
plane obtained by Monte Carlo simulation (circles) and by 
mean field theory (dashed curves). The Higgs-Coulomb and 
confinement- Coulomb phase transitions are of second order, 
whereas the Higgs- confinement transition is of first-order in 
the region near the tricritical point, and becomes second-order 
for higher ci. For example, (c2,ci) = (0.96,1.27) is of first 
order and (0.87,1.40) is of second order. 



in the time direction ao suitably, we obtain Eq.([Tj) with 
ci a constant independent of the exchange coupling J of 
the AF Heisenberg model, and C2 = 0[3]. However, via 
the renormalization effect of high-momentum modes of 
spinous not only the ci-term is renormalized but also 
the C2-term is generated. 

When holes are doped into an AF magnet like the high- 
Tc cuprate at half-filling, the system is described by a 
canonical model like the t-J model. The t-J model can be 
studied by the slave- fermion- CP ^ representation [1]. One 
can see that as a result of the hole doping in the short- 
range AF background, the parameter ci is renormalized. 
If the hopping of holons is ignored, ci {1 — S)ci where 
5 is the hole concentration. Furthermore, the short-range 
AF order generates attractive force between holes sitting 
on nearest-neighbor sites, and it indicates the appearance 
of the superconducting phase. When the hole-pair field 
Mx^i (i = 1,2,3 : spatial direction index) condenses, the 
C2-plaquette term in Eq.([T|) is generated by the hopping 
of quasiparticles (i.e., the gapped holons) [1]. 

The phase diagram of the model ([T]) in the C2 — ci 
plane has been determined by calculating the "internal 
energy" per site E = {S)/V and the "specific heat" per 
site C = {{S — EY)/V by means of Monte Carlo simula- 
tions for a lattice of size V = with the periodic bound- 
ary condition in the previous paper [4]. We show the re- 
sult in FiglH In FiglU we also show the result obtained 
by the mean-field approximation. There are three phases 
in the C2 — ci plane. We expect and verify that these three 
phases are Higgs, confinement and Coulomb phases in the 
gauge-theory terminology, respectively. The Higgs phase 
is nothing but the Neel state of the AF magnets. There 
deconfined gapless spinous are spin- wave magnons. The 
confinement phase is a spin-liquid phase in which the low- 
energy excitations are spin-triplet hound states of Zx- On 
the other hand, the Coulomb phase corresponds to a de- 
confined spin liquid in which the low-energy excitations 
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FIG. 2: Spin-spin correlation functions Cs(t) for 12^ lattice, 
(a) C2 = 2.0; (b) ci = 1.5. The magnetization is finite in the 
Higgs phase. 



are s = ^ spinous Zx with a gap. We shall see that the 
gauge-boson mass vanishes in the Coulomb phase while 
it remains finite in the other two phases. 

In Ref.[5] we have studied the phase structure of the 
CP^+U(1) model ^ not on the 4D lattice but on the 
3D lattice. We found that the Coulomb phase is missing. 
The gauge-boson mass takes its minimum value just on 
the critical line separating the Higgs and confinement 
phases. If we include A^f-fold CP^ variables, the gauge- 
boson mass on the critical line vanishes for A^f > 14 [y]. 
The original uniform AF Heisenberg model is in the Neel 
state, which corresponds to the Higgs phase [1]. 

In Fig|2l we show the result of our calculations of the 
0(3) spin-spin correlation functions, Cs(r) = {Sx+r'Sx)- 
From Fig|2] we see that the spontaneous magnetization, 
[Cs(rmax)]^/^ (^max = i^), is nouvauishing only in the 
Higgs phase as expected. 

In Figl3]we show the mass gap Ms of the spin excita- 
tions that is obtained from Cs(r) — Cs(rmax)- Fig |3] shows 
that Ms is vanishing in the Higgs phase as expected. 

Next, we study the dynamics of gauge bosons by mea- 
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FIG. 3: Mass Ms of spin excitations for C2 = 2.0. P.T. indi- 
cates the location of transition point in Fig|T] Ms vanishes in 
the Neel state (Higgs phase). 
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FIG. 4: Gauge-boson mass Mq of (gj). (a) ci = 1.5 (b) ci = 
0.5 (c) C2 = 2.0. Mg vanishes only in the Coulomb phase. 
(Negative value of Mg is a finite-size effect. See Ref.[6].) 



suring the gauge-boson mass Mq, which is calculated 
from the correlation function of gauge flux, 

where e^^xr^ is the totally antisymmetric tensor. To esti- 
mate Mq precisely, we introduced the Fourier-transform 
O and measured its correlation, 



)),(3) 



with setting px = Prj = l/L where L is the system size. 
We fit the data in the following form^]; 



Dciy^,, Vu) (X exp [-^ pI ^ yi) . (4) 

In FigHl we show Mq. We can see that Mq is vanishing 
in the Coulomb phase as it should be, whereas it is finite 
in the other phases. 
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FIG. 5: AF Heisenberg model with dimer and ring couplings. 
The double-line bonds have larger exchange than that of 
the single- line bond. The wavy lines represent ring coupling 
among four spins on each square of even-sites. 
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FIG. 6: Specific heat of the 3D plaquette model (|5]) for ci = 
0.5 as a function of S2- There is no system-size dependence 
in C, indicating a crossover near the round peak. 



The above measurements confirm the phase structure 
of the model ([1]) given in Fig{TJ From this result, we can 
discuss how the magnetic properties of the doped AF 
magnets with strong three-dimensionality changes under 
doping at zero temperature. Undoped and hghtly doped 
uniform AF magnets in 3D has the Neel order and is 
in the Higgs phase. As 5 is increases, the system loses 
the Neel order and enters the confinement or Coulomb 
phase depending on the mobility of doped holes. In the 
framework of the t-J model, this problem is under study 
and the result will be published in near future [8]. 

In the rest of this paper, we shall consider the CP^ 
gauge model with the 3D spatial gauge-plaquette term. 
The action S' of this 3D plaquette model is given by 



S' = 



2 ^ ^ (^x-\-n^aUx^i^xa H~ H.C.^ 

2 

3 

(UxjUx+j,iUx+i,jUxi + H.c.^(5) 



£2 

2 



x,ij=l{i<j) 

In contrast to the C2 term of the 4D plaquette model ([T]), 
the S2-term in Eq.([5|) contains only the spatial-plaquette 
terms and no space- time plaquette terms. The 3D pla- 
quette model ([5|) is an effective low-energy model of 
nonuniform AF Heisenberg model with dimer-like cou- 
plings and ring-exchange couplings. The parameter ci 
decreases as the dimer coupling increases [9] and the 52 
term corresponds to ring coupling like (5'^ • Sj){Sk • Si) 
where (i, j, /c, /) denote even sites forming a square in even 
sublattice. (See Fig|5l) Thus the bare value of 52 is finite 
in contrast to the 4D plaquette model ([1]) for the sym- 
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FIG. 7: Internal energy of the 3D plaquette model (|5]) for 
S2 = 5.0. It shows a hysteresis curve. 
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FIG. 8: Gauge boson mass Mg in the model (|5]). (a) ci — 0.5; 
(b) S2 = 5.013. 



metric Heisenberg model where the bare value of C2 = 0. 

It is interesting to see if the Coulomb phase (decon- 
fined spin-liquid phase), which exists in the model ([T]), 
survives in the present model. Therefore we investigated 
the phase structure of the model ([5]) first by calculating 
the internal energy and the specific heat. Some of our 
results are shown in FigslH and [71 From these calcula- 
tions, we conclude that the model S' of ([5]) has similar 
phase structure to that of 5 of ([T|). But the second-order 



confinement- Coulomb phase transition in the model S 
becomes a crossover in S' as Fig l6] shows, i.e., the spe- 
cific heat has a rather smooth peak that has no system- 
size dependence. We think that this crossover is similar 
to that in the 3D Abelian gauge-Higgs model. Further- 
more, the phase transitions at ci :^ 0.7 in the region of 
82 ^ 1.5 are of first order as indicated in Fig 171 

In order to verify the above conclusion, we calculated 
Mg, which is shown in FiglSl It is obvious that along 
with ci = 0.5 Mg decreases smoothly as 52 is increased 
crossing the crossover line at 52 — 1.5. This behavior 
is very close to that in the Abelian gauge-Higgs system 
observed in Ref.[6]. On the other hand, along the line 
S2 = 5.0 Mg exhibits a sharp discontinuity. Its behavior 
verifies the first-order phase transition at ci O.Tpioj. 

Then we draw the following conclusions; 
1. In contrast to the 4D plaquette model the 3D pla- 
quette model does not support the Coulomb phase, 

1. e., the deconfined spin-liquid phase is not realized in 
the nonuniform AF Heisenberg model [11, 12]. 

2. The Neel-VBS phase transition caused by the ring- 
exchange coupling is of first-order. Recent numerical 
study on some related AF Heisenberg models gives simi- 
lar conclusion [l2|. 

In summary, we studied the 4D CP^ model coupled 
with 4D or 3D U(l) gauge field and revealed its phase 
structure. We discussed some physical implications for 
various 3D AF Heisenberg magnets, i.e., with and with- 
out doping and with uniform or nonuniform coupling in 
gauge-theoretical viewpoints. 
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